HD gas analysis with Gas Chromatography and Quadrupole Mass Spectrometer 
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Abstract 

(N ■ 

A gas analyzer system has been developed to analyze Hydrogen-Deuteride (HD) gas for producing frozen-spin polarized HD 
^ targets, which are used for hadron photoproduction experiments at SPring-8. Small amounts of ortho-H2 and para-D2 gas mixtures 
H-} (~0.01%) in the purified HD gas are a key to realize a frozen-spin polarized target. However, there was an intrinsic difficulty 
qq to measure these small mixtures in the HD gas with a quadrupole mass spectrometer (QMS) because D + and [H2D] + produced 
^vq from the ionization of HD molecules were misidentified as H2 and D2 molecules, respectively, and became backgrounds for the 
measurement of the H2 and D2 concentrations. In addition, the ortho-H2 and para-D2 are not distinguished from the para-H2 and 
' ^ ' ortho-D2, respectively, with the QMS. In order to obtain reliable concentrations of these gas mixtures in the HD gas, we produced 
dJ a new gas analyzer system combining two independent measurements with the gas chromatography and the QMS. Helium or neon 
1 ( gas was used as a carrier gas for the gas chromatography which was cooled at ~ 1 10 K. The para-H2, ortho-H2, HD, and D2 are 
O separated using the retention time of the gas chromatography and the mass/charge ratio. Although the para-D2 is not separated from 
^ the ortho-D2, the total amount of the D2 is measured without the [H2D] + background. The ortho-H2 concentration is also measured 

1 ^ 1 separately from the D + background. It is found that the new gas analyzer system can measure small concentrations of ~0.01% for 

the otho-H2 and D2 with good S/N ratios. 
y—i ■ 
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IO !• Introduction 

Ln ; 

__, Double polarization measurements for <p and K meson photo 
productions with a polarized target and a polarized photon beam 

7 ~ I are a sensitive means to investigate hadron structures, such as 
I ian ss-quark content of nucleons |l}] and reaction mechanisms 
!>" jSSBl- In order to measure such "complete polarization ob- 

"k^j servables", an essential key technology is to prepare a polarized 
target with a long relaxation time. We have been developing 
C\3 a frozen-spin polarized Hydrogen-Deuterium (HD) target on 
the basis of a new technologies in cryogenics, superconducting 
magnet, and HD gas purification. These developments would 
lead us to the next upgrades of the LEPS experiments [51]. An 
idea for making the frozen-spin polarized HD target came from 
the work in Refs. It is suggested that the HD molecule 

can be used as a frozen spin polarized target [8]. Many ef- 
forts have been devoted for a long time to realize the polarized 
HD target at Syracuse H, [H, BNL id [H 0, and ORSAY 
EH HE Gj, OH. Recently, the HD target has been firstly 
used for the actual experiment at LEGS 01911 . and will be used 
at JLab [20] and at SPring-8 J2l| in near future. 
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In order to prepare the frozen-spin polarized HD target, 
Honig [8] has applied an innovative idea that the HD polariza- 
tion gradually grows up by the spin-flip process between HD 
molecules and a small amount of ortho-H2 (0-H2) with spin 
J=l, whose aspect was many year ago suggested by Motizuki 
et al. I22I E3I1 . Meanwhile after a long period of cooling, 0-H2 
molecules are converted to para-H2 (P-H2) with spin J=0. If all 
the 0-H2 molecules are converted to the P-H2 in the HD target, 
the relaxation time of hydrogen polarization becomes very long 
even at a temperature higher than 4 K. After an aging process 
for 2-3 months, the long relaxation time is expected with a hy- 
drogen polarization higher than 84%. The polarization degree 
is measured by using the NMR method l24ll . In the conversion 
process from 0-H2 to P-H2, an important parameter is a quan- 
tity of remaining 0-H2 with J=l, which acts as a mediator for 
spin depolarization when the temperature of the solid HD target 
rises up to more than 4 K: the lower the impurity of 0-H2, the 
longer the relaxation time of hydrogen. Thus, it is essential to 
control the admixture of the 0-H2 in the highly purified HD gas. 

When the HD gas is kept in the gas container for a long time, 
the quality of the HD gas is deteriorated due to the dissocia- 
tion process, 2HD < — > H2 + D2 and the H2 and D2 gasses are 
naturally yielded. For the HD purification, we installed a distil- 
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lation still to make the pure HD gas with a purification level of 
99.99%. The next indispensable problem to be solved is how to 
measure the small admixtures of 0-H2, P-H2, and D2 in the HD 
gas. The conventional method is to use the quadrupole mass 
spectrometer (QMS) for gas analysis. This instrument makes 
use of a principle that the analyzed gas is ionized at an ion 
source by electron bombardment and ions are mass-separated 
according to the mass/charge ratio (u/e). When the huge HD 
gas is entered in the ion source of the QMS, small amounts of 
H + and D + are produced from HD molecules. Other ions of H^ , 
[H2D] + , [HD2] + , and D3 are also produced by recombination. 
This means that D + and [H2D] + are misidentified as HJ and DJ , 
respectively. The QMS cannot distinguish the molecules and 
fragments with a same mass/charge ratio. We call this prob- 
lem as "fragmentation problem". The fragments produced by 
ionization are listed in Table Q] 



the gas is very long. On the other hand, since the hydrogen and 
deuterium are light and small particles, the retention times of 
0-H2, P-H2, HD and D2 are as short as a few second at the room 
temperature. Because of this, the zeolite column is cooled down 
to about 1 10 K to attain a reasonably long retention time for the 
hydrogen and deuterium gases. Actually, the column is cooled 
down in a constant temperature dewar by controlling the vapor 
flow from liquid N2 (LN2). The capillary column is installed 
in the dewar with the vacuum insulation layer. LN2 is stored at 
the center layer, and the cold vapor gas from LN2 flows from 
the inner layer to the outer layer. We use the He or Ne gas as a 
carrier gas. The flow rate of the gas is controlled by using a flow 
controller. The sample which is out of the column is detected 
by the QMS. This QMS is a product of MKS Instruments JH 
and is specially tuned-up to increase the sensitivity for detecting 
light molecules such as H2, HD, D2, and T2. 



Table 1: The isomers in the HD gas and expected fragments. 



u/e 1 


2 


3 


4 


5 6 


Molecules 


H 2 + 


[HD] + 


D+ 




Fragments H + 


D + 


H| 


[H 2 D] + 


[HD 2 ] + D+ 



To overcome the fragmentation problem, we designed a sys- 
tem combined with the gas chromatography and the QMS for 
the HD gas analysis with a high dynamic range of more than 
10 4 . The gas chromatography separates the isomers in the HD 
gas in terms of time. The QMS separates the isomers in terms of 
mass/change. By combining two different analysis systems, it is 
possible for us to analyze the H2 and D2 concentrations in the 
HD gas, precisely. The James Madison University group has 
successfully employed a gas chromatography for distinguish- 
ing 0-H2 and P-H2 [25]. However, a small D2 admixture with a 
level of 0.01% cannot be clearly observed because an HD long 
tail conceals a tiny D2 peak in the gas chromatography spec- 
trum. 

We have developed a new gas analyzer system by combining 
the gas chromatography and the quadrupole mass spectrometer 
for removing these disadvantages. Precise measurements of the 
0-H2 and D2 concentrations in the HD gas will be used to effec- 
tively polarize the HD target and to obtain the long relaxation 
time. 

2. System overview 

Fig. [TJ shows a schematic diagram of a new gas analyzer 
system. The HD gas with a volume of 500 fil is infused by a 
gas sampler into the column. The sample gases are pushed in 
the column together with a carrier gas. The isomers of 0-H2, 
P-H2, HD, and D2 in the HD gas are separated by using the 
difference of adsorbent action to zeolites which are affixed on 
the inner wall surface of the column. The separated isomers of 
the HD gas are analyzed by the quadrupole mass spectrometer 
as a function of the mass/charge ratio. In the usual method of 
a gas chromatography, the temperature of a thin zeolite column 
is increased during the analysis because the retention time of 
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Figure 1: Schematic diagram of the gas chromatograph and the quadrupole 
mass spectrometer for the HD gas analysis. The HD gas is injected with a 
carrier gas from the entrance of the column via gas sampler. QMS: quadrupole 
mass spectrometer. TMP: Turbo molecular pump. 



2.1. Distillator 

To make a polarized HD target with a long relaxation time 
of more than one month, we need to prepare a highly purified 
HD gas. However, if we use highly purified HD gas, the aging 
time for polarizing the HD becomes long. This is an experi- 
mental dilemma. One of methods to overcome this dilemma is 
to dope an appropriate amount of 0-H2 and P-D2 impurities in 
the purified HD gas. The appropriate amount is approximately 
an order of 0.01%. The purity of the commercially available 
HD gas is about 96%. Contaminations are mostly H2 (~ 2%) 
and D 2 (-2%). 

The HD gas is necessary to be purified up to > 99.99% for 
optimizing the amount of impurities. The gas distillation sys- 
tem used to purify the HD gas at RCNP is shown in Fig. [2] The 
commercial gas is fed to the pot inside the ditillator. There are 
stainless cells called "Helipack" inside the pot. Temperature 
gradient is realized by cooling the top of the pot and heating the 
bottom by the thermal radiator. Heat exchange between gas and 
liquid takes place on the cells. A gas with a low boiling tem- 
perature is extracted from the top of the distillator. By using the 
difference of the boiling temperatures of HD (16.6 K), H2 (14.0 



2 



K) and D2 (18.7 K), we separate the HD gas from the others. 
In the initial trial, the concentrations of HD, H2, and D2 were 
measured with the QMS. Although the H2 and D2 concentra- 
tions were decreased less than 0. 1 %, it was not able to measure 
the concentrations precisely because of the fragmentation prob- 
lem. 




Figure 2: (a) The design view of the distillator. The H2 gas is extracted from 
the top of the distilllation pot. (b) The photo picture of stainless steel packs 
which are filled into the pot in the distillator. Each pack has a rolled shape like 
a coil and twisted for efficient heat conduction between gas and liquid. 



can't pass through the quadrupole rods because of the resonance 
mismatch. The separated ions are detected at the Faraday cup as 
electric current from a multiplier. All isomers in the HD gas can 
be measured by scanning the applied frequency, independently. 
We operated the QMS with a pressure of about 10 9 -10 8 bar. 
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Figure 3: Structure of the quadrupole mass spectrometer. Constant (U) and al- 
ternative (V cos(wt)) voltages are applied to four cylindrical electrodes. The in- 
put HD gas is ionized in the ionization chamber by electron bombardment. Ions 
are injected along the central symmetry axis of the four cylindrical electrode- 
rods. An ion satisfying a resonance condition with a frequency to makes a 
helical movement coupled to a displacement movement towards the Faraday 
cup, and is finally detected with a multiplier. 



2.3. Gas chromatography(GC) 

We use the Molsieve 5A PLOT (porous layer open tubular) 
type column. The capillary column is made of fused silica and 
coated by polyimide. The size of the capillary column is sum- 
marized in Table 2. 



We operated the distillator for one week to obtain pure HD 
gas. The commercial HD gas with an amount of 5.3 mol was fed 
to the pot. A gas was extracted from the distillator with a flow 
rate of 2.0 ml/min to the gas storage tanks made of stainless 
steel. The H2 concentration was very high in the beginning 
and the HD purity gradually increased in a few days. After 
extracting about 1 mol gas, the gas analysis started. 

2.2. Quadrupole mass spectrometerf QMS) 

A quadrupole mass spectrometer (QMS) consists of an ion 
source, an extraction plate, four cylindrical electrodes, and a 
Faraday cup. A high mass resolution and compactness are its 
good features. The present QMS (made by MKS Instruments, 
Microvision Plus |26J]) is equipped with a multiplier to mea- 
sure a partial pressure down to 10~ 17 bar. The QMS has an abil- 
ity to well distinguish the masses of He (amu=4.0026) and D2 
(amu=4.0282), because its performance is customized to mea- 
sure low mass molecules such as H2, D2, and T2 in the mass 
range of 1-6 Q. 

Constant and alternative voltages are applied to four cylin- 
drical electrodes (quadrupole rods) as shown in Fig. [3] The 
sample gas is ionized in an ionization chamber by an elec- 
tron impact method and then accelerated by the extraction plate 
along to the symmetry axis of the quadrupole rods. Ionized 
ions move towards the Faraday cup with a spiral motion cou- 
pled to a displacement motion. Ionized ions, which satisfy the 
resonance condition of the Mathieu's differential equation, pass 
the quadrupole rods and reach the Faraday cup. The other ions 



Table 2: Specification of the fused silica column. Molsieve 5A PLOT type 
column 1 28] was used. OD and ID stand for the outer and inner diameters of 
the column. 

OD ID Length Zeolite thickness 

0.70 mm 0.50 mm 50 m 0.050 mm 



Fig.|4]shows an expected spectrum of gas chromatography as 
a function of retention time. The peak in the gas chromatogram 
data was analyzed by fitting with the exponential-Gaussian hy- 
brid function B29I1 as follows: 



fit) 



Hexp 
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where H is the peak height, cr f is the standard deviation of the 
Gaussian, t is the time constant of the exponential, and t R is the 
time of the peak maximum. 

We define the separation degree R between two peaks as, 

(fz-h) 



R 



(2) 



(0"i + cr 2 ) 

where tj and t2 are their retention time, cr\ and 1x2 are calculated 

as, 



0"1 



-4 



^1 + A, 0-2 



-4 



(3) 



where <r g \ and cr g2 are the standard deviation for the peak 1 and 
the peak 2, respectively, and t\ and r 2 are the time constant for 
the peak 1 and peak 2, respectively. 

Longer column length and longer retention time are consid- 
ered to give large separation R. But very long retention time is 
not appropriate for the present work. To measure the HD gas 
from the distillator efficiently, we have estimated the best R by 
changing various parameters such as the temperature of the col- 
umn and the flow rate of the He carrier gas. 




Figure 4: Expected peaks in the gas chromatography. t\ and ?2 are the retention 
times, ri and T2 are the time constant of the peak. The separation degree R is 
defined by Equationf2] Two peaks are separated completely with R > 1.5. 



3. Experimental results and analysis 

3.1. Effect of temperature and gas flow 

For precise measurement of the concentration of isomers in 
the HD gas, R is required to be more than 1.5. To determine the 
optimum temperature and the flow rate, the separation degree 
R between HD and D 2 was measured by changing various pa- 
rameters of the temperature and the flow rate as shown in Table 
[3] Lower temperature (105K) with a flow rate of 10ml results 
in the best separation. 

3.2. Measurement of combination with GC and QMS 

We prepared a sample gas by mixing appropriate amounts of 
H 2 , HD, and D 2 gases. The sample was analyzed by the gas an- 
alyzer system. The neon gas with a flow rate of 35 ml/min was 
used to carry the sample through the column cooled at 105 K. 
The p-H 2 and o-H 2 were detected at 8 and 9 minutes, respec- 
tively, after the injection of the sample as shown in Fig. EJa). 
The HD and D 2 gases were detected at 1 1 and 14 min, respec- 
tively, as shown in Fig. Qb) and (c). Using both the GC and 
the QMS, we could observe the p-H 2 , o-H 2 , HD, and D 2 gases 
separately by measuring the retention time in the GC and by 
determining the mass/charge ratios in the QMS as shown in 
Fig. [6] When the purity of the sample HD gas becomes very 
high (~99.9%), backgrounds disturb clear obsevation of the H 2 
and D 2 peaks. This problem is discussed in Sec. 13.41 




Elapsed Time (min) 



Figure 5: Gas chromatograms for the gases with the mass/charge ratios of (a) 
2 (b) 3 and (c) 4 . The vertical axis is the partial pressure for each gas and the 
horizontal axis is the elapsed time after the injection of the gas sample. 




Figure 6: Three-dimensional plot of the elapsed time of the GC vs the 
mass/charge ratio measured by the QMS. The z axis is the partial pressure of 
each gas. 
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Table 3: The separation degree R between HD and D2. Errors are the sum of the statistical and systematic errors. 



Flow rate 


1 ml 


3 ml 


5 ml 


10 ml 


20 ml 


30 ml 


105 K 
110K 
120 K 
130 K 


2.610+0.002 
2.369±0.002 


4.881 ±0.005 
4.383+0.018 
2.416+0.004 


5.754+0.009 
5.346±0.032 
7.425+0.054 
3.946+0.087 


7.988+0.042 
7.646+0.017 
7.702+0.061 
5.253+0.027 


4.568+0.007 
6.107+0.013 
6.800+0.043 
3.748+0.062 


5.346+0.008 
6.685+0.026 
6.862+0.055 
3.802+0.016 



3.3. Transition between p-H 2 and o-H 2 

There are mysterious events between the P-H2 and 0-H2 
peaks in the gas chromatograms as shown in Fig. |5ja). These 
events are considered to be produced by the transition from p- 
H2 to 0-H2 or 0-H2 to P-H2 in the column at low temperature. If 
the P-H2 or 0-H2 molecule interacts with magnetized materials, 
the transition can be induced although the materials have not 
been clearly specified yet. Such events between the P-H2 and 
0-H2 peaks are observed in another experiment ll25ll . 

We fit the gas chromatogram data as shown in Fig.|7]with the 
function consisting of two exponential-Gaussian hybrid func- 
tions written in Eq. Q] a function (Eq. |4]i for reproducing the 
transition events written below, and a constant for background 
events. 
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where C is a normalization factor. Hi and H 2 are the p-F^ and 
0-H2 peak heights, cr gl and cr g2 are the standard deviation of the 
Gaussian, t\ and t 2 are the time constants, and 1 1 and t 2 are the 
time of the peak maximum, respectively. 

The errors for the output of the QMS are estimated from the 
measurement fluctuation at a stable pressure. The reduced x 1 
for the fit is 11.3. The tail of the 0-H2 peak is not well repro- 
duced, which makes the reduced^ 2 large. The concentrations 
of P-H2, 0-H2, and transition events are obtained as shown in 
Table|4] The errors of the parameters in the fit and the deviation 
between the fit and the data at the 0-H2 peak tail are considered 
as errors of the concentrations. 

3.4. Analysis for pure HD 

The best separation of the peaks was obtained when the he- 
lium carrier gas with a flow rate of 10 ml/min was used at 105 
K as listed in Table [3] When the purity of a sample HD gas 
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Figure 7: The fitting to the data with 2 u/e. Peaks of 0-H2 and P-H2 are fit by 
Eqn.fT] Transition region is fit by Eqn.|4] 



Table 4: Concentrations of p-Eb, 0-H2, and transition events obtained by the fit 
to the data. 





para-H2 


para-H2 » ortho-H2 ortho-H2 


Concentration 


30+1% 


6±1% 63±2% 



becomes very high, the experimental condition needs to be op- 
timized in order to reduce backgrounds. A pure sample HD 
gas was analyzed by using a helium carrier gas with a flow rate 
of 1.0 ml/min at 125 K. As shown in Fig. IS a), the HD peak 
is dominantly observed at 16.5 min in 3 u/e. Another peak is 
also observed at the same position in 2 u/e. This peak is due 
to D + produced from the ionization of HD. Without the GC, 
the D + component could not be easily separated from real H2 
signals and the measurement of the H2 concentration was un- 
certain. The 0-H2 peak is observed at 15.5 min separately from 
the P-H2 peak at 14.1 min. The concentrations of P-H2 and 
o-H 2 were obtained as 0.005+0.001(stat.)+0.001(syst.)% and 
0.010+0.001(stat.)+0.002(syst.)%, respectively. The H 2 events 
between the P-H2 and 0-H2 peaks and the tail of the peaks were 
assumed to be the P-H2 or 0-H2 with a ratio of 5 (P-H2) to 10 
(0-H2). Present background level in 2 u/e gas chromatogram 
enables us to measure the 0-H2 concentration with a precision 
better than 0.01%. 

The linearity of the output of the QMS and the detection ef- 
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ficiency of the GC were taken into account as the systematic 
errors. We checked the linearity in the partial pressure region 
of 10~ 16 -1(T 12 bar. The relation between the amount of sam- 
ple gas and the partial pressure measured by the spectrometer 
were well fitted with a linear function. The deviation from the 
linear function was smaller than 12% of the measured partial 
pressure. The detection efficiency was measured by storing the 
sample gas and the carrier gas in a tank after passing through 
the gas chromatography and the QMS. The uncertainty of the 
detection efficiency was found to be about 17%. 

The pure sample HD gas was analyzed by using a neon car- 
rier gas with a flow rate of 1.0 ml/min at 125 K. As shown in 
Fig.[8jb), the HD peak is dominantly observed at 16.0 min in 
the gas chromatogram with 3 u/e. At the same location of the 
HD peak, another peak is observed in 4 u/e. This peak is in- 
ferred to appear due to the [H2D] + molecule produced by the 
ionization of HD, which was confirmed by the result that the 
[HD2] + was also observed at the same location in the gas chro- 
matogram with 5 u/e (not shown). The D2 peak is clearly ob- 
served at 19.8 min on top of the long tail of the [H2D] + peak. 
The concentration of D2 in the sample HD gas was obtained 
as 0.043+0.001(stat.)+0.009(syst.)%. Judging from the fluctu- 
ation of the background around the D2 peak, the D2 concen- 
tration of 0.01% can be detected by using the present new gas 
analyzer system. 



Table 5: The measured concentrations of p-tb, 0-H2, HD, and D2 in the dis- 
tilled HD gas 
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Isomers 


Concentration 


P-H 2 


0.005 ±0.001(stat.) +0.001(syst.) % 


o-H 2 


0.010 ±0.001(stat.) +0.002(syst.) % 


HD 


99.942 ±0.002(stat.) +0.009(syst.) % 


D 2 


0.043 ±0.001(stat.) +0.009(syst.) % 



Figure 9: (a) The gas chromatogram with 4 u/e for a distilled HD gas measured 
by using the neon carrier gas. (b) The sum of the gas chromatograms with 3 
u/e and 4 u/e for a distilled HD gas measured by using the neon carrier gas. 
The vertical axis is the partial pressure of each gas. The horizontal axis is the 
elapsed time. 



In Fig. |9] a D2 peak in gas chromatogram with 4 u/e, mea- 
sured by using the gas chromatography and QMS, is compared 
with another small D2 peak in the gas chromatogram, measured 
by using the gas chromatography only. The D2 peak is clearly 
observed with a good S/N ratio in Fig.|9ja). However, it is very 
difficult to find the small D2 peak because of the large back- 
ground from the tail of the HD peak in Fig. |9jb). The S/N ratio 
for the D2 peak was improved by a factor of 10 in the measure- 
ments by using both the GC and the QMS. 



concentrations of the P-H2, 0-H2, and D2 in the HD gas is de- 
voted at JLab [30]. The JLab group is trying to analyze the 
HD gas by using the Raman scattering of laser light. The ac- 
curacy of measuring the concentrations will be improved by 
introducing the Raman spectroscopy although the setup is aw- 
fully complicate. It should be noted that the developments of 
gas analyzing techniques reported in the present work will play 
an important role in producing the polarized HD target under a 
well defined HD gas condition. 



4. Summary 

We have developed a new HD gas analyzer system by com- 
bining the gas chromatography and the quadrupole mass spec- 
trometer for producing the polarized HD target to be used 
for the hadron photoproduction experiments at SPring-8. The 
new system enabled us to observe P-H2, 0-H2, HD, and D2 
separately. We succeeded in measuring small concentrations 
(-0.01%) of p-H 2 , o-H 2 , and D 2 in the distilled HD gas with 
good S/N ratios. Recently, another effort to measure the small 
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Figure 8: (a) Gas chromatograms for a distilled HD gas measured by using the helium carrier gas. The gas chromatograms with 2 u/e (solid curve) and 3 u/e (dotted 
curve) are plotted, (b) Gas chromatogram for the distilled HD gas measured by using the neon carrier gas. The gas chromatograms with 3 u/e (dotted curve) and 4 
u/e (solid curve) are plotted. The vertical axis is the partial pressure of each gas. The horizontal axis is the elapsed time. 
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